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We investigate the electrodynamic structure of a pulsar outer-magnetospheric particle accelerator and the resul- 
tant gamma-ray emission. By considering the condition for the accelerator to be self-sustained, we derive how 
the trans-magnetic-field thickness of the accelerator evolves with the pulsar age. It is found that the thickness 
is small but increases steadily if the neutron-star envelope is contaminated by sufficient light elements. For 
such a light element envelope, the gamma-ray luminosity of the accelerator is kept approximately constant as 
a function of age in the initial ten thousand years, forming the lower bound of the observed distribution of the 
gamma-ray luminosity of rotation-powered pulsars. If the envelope consists of only heavy elements, on the other 
hand, the thickness is greater but increases less rapidly than what a light element envelope has. For such a 
heavy element envelope, the gamma-ray luminosity decreases relatively rapidly, forming the upper bound of the 
observed distribution. The gamma-ray luminosity of a general pulsar resides between these two extreme cases, 
reflecting the envelope composition and the magnetic inclination angle with respect to the rotation axis. The 
cutoff energy of the primary curvature emission is regulated below several GeV even for young pulsars, because 
the gap thickness, and hence the acceleration electric field is suppressed by the polarization of the produced 
pairs. 



1. Introduction 

The Large Area Telescope aboa r d Fer mi Gamma- 
ray Space Telescope ( Atwood et al.l [2009j ] ^ has proved 
remarkably successful at discovering rotation-powered 
pulsars emitting photons above 0.1 GeV. Thanks to 
its superb sensitivity, the number of gamma-ray pul- 
sars has increase d from six in Com pton Gamma Ray 
Observat ory era (iThompsonl |2004j ) to more than one 
hundred (|Nolanl |2012j ). Plotting their best estimate 
of the gamma-ray luminosity, L 
down luminosity, L ; 
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4tt 



7 , against the spin- 
2 JPP- 3 , lAbdo et alJ 



[201C| found the important relation that L 1 is approx- 
imately proportional to L sp i n ' 5 (with a large scat- 
ter), where I refers to the neutron-star (NS) mo- 
ment of inertia, P the NS rotational period, and P 
its temporal derivative. However, it is unclear why 
this relationship arises, in spite of its potential im- 
portance to discriminat e pulsar emission models such 
as the polar-cap model dDaughertv fc Harding] [l982] , 
iDermerl fl994ll. lHarding et all 1197811 the ou t er-gai 
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model (I Chian g fc Romanil 

T996t . ITakata et all 



2011| . IZhang fc Chend Il997l 



1992 1. iHirotanil 1200 



2004} , 
the 
2009]) 



Wang et all 



pair-starved 
(see also 



polar-cap mod e l (IVenter et all _ 
lYuki fc Shibatal [2012j for the possible co-existence 
of such m odels), and the emissio n model from the 
wind z on e (lAharonian et all [2012 j , iBai fc Spitkovskil 
(2010allbl |. |Petril [201 l| l 

Recent gamma-ray observations suggest that the 
pulsed gamma-rays arc emitted from the higher al- 
titudes of a pulsar magnetos phere. Thi s is be cause 
the observed light curves (|Abdo et all |2010j ) fa- 
vor fan-like emission geometry, which scan over a 
large fraction of the celestial sphere, and because 
the Crab pulsar shows pulse d photo ns near and 
above 100 GeV (lAleksic et al.l [2011 alibi ]. lAliu et all 



201 1| ). which rules out an emission from the lower 
altitudes, where strong magnetic absorption takes 
place for 7-rays above 10 GeV. Consequently, higher- 
altitude emiss i on mod els such as the outer-gap model 
( Cheng et al. Il986allbl0. the high-a ltitude slot-gap 



model ( Muslimov fc Harding! 1200411 or th e pair- 
starved polar-cap model ( Venter et al.l [20091 ] ). gath- 
ered attention. It is noteworthy that the outer-gap 
model is presently the only higher-altitude emission 
model that i s solvable from t he basic equations self- 
consistently (|Hirotanil [2011aj l In the present paper, 
therefore, we focus on the outer-gap model and derive 
the observed relationship L 7 oc L spnl 5 both analyti- 
cally and numerically. 



We schematically depict the pulsar outer- 
magnetosphe ric accelerator (i.e., the outer gap) 
in figure 1 of lHirotanil [20131 ] As the NS rotates, there 
appears the light cylinder, within which plasmas can 
co-rotate with the magnetosphere. The magnetic field 
lines that become tangential to the light cylinder at 
the light cylinder radius, tolc = cP/2n, are called the 
last-open magnetic field lines, where c refers to the 
speed of light. Pairs are produced via photon-photon 
pair production mostly near the null-charge surface 
and quickly polarized by the magnetic-field-aligned 
electric field, , in the gap. In this paper, we 
assume that the rotation and magnetic axes reside 
in the same hemisphere to obtain E\\ > 0, which 
accelerates positrons (e + 's) outwards while electrons 
(e~'s) inwards. These ultra-relativistic particles have 
Lorentz factors, 7 ~ 10 7 ' 5 , to emit photons efficiently 
by the curvature process. 
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2. Analytical examination of outer-gap 
luminosity 

In this section, we analytically derive how the 
gamma-ray luminosity of an outer gap evolves with 
time. In the outer magnetosphere, only the dipole 
component remains in the magnetic field; thus, the 
inhomogeneous part of the Maxwell equation (i.e., 
the Poisson equation for the electro-static potential) 
gives the magnetic-field-aligned electric field [Hirotanil 
20081. 
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(1) 



where fi denotes the NS magnetic dipole moment, 
and h m the trans-magnetic-field thickness of the gap. 
Since the Poisson equation is a second-order differen- 
tial equation, E\t is proportional to h m 2 . Electrons 
(e~'s) and positrons (e + 's) are created via photon- 
photon (and sometimes via magnetic) pair produc- 
tion, being subsequently polarized by E\\ and acceler- 
ated in the opposite directions, to finally attain the 
terminal Lorentz factor 



7 



2e 



1/4 



where p c refers to the radius of curvature of parti- 
cle's motion in the three-dimensional magnetosphere, 
e the charge on the positron. Photons are radiated by 
such ultra-relativistic e^'s via curvature process with 
characteristic energy, 



hv c = -Tic — , 
2 p c 



(3) 



where h denotes the Planck constant, fi = h/2n. Once 
h m is obtained, we can readily compute the 7-ray lu- 
minosity of cur vature radiation from an outer gap by 
[Hirotanil [2008j 



L 7 «2.36(^) pcak x47rd 2 /fi 
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l.A6fnti m — r— . 
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where /n, which has been conventionally assumed 
to be appro ximately unity, refers to the flux correc- 
tion factor IRomani. R. fc Watterll2010j . and = 
2ir/P the rotation angular frequency of the NS. Here, 
it is assumed that the current density flowing in 
the gap is comparable to the Goldreich- Julain value 
Goldreich. fc Julian! fl969lj . The last factor, p 2 fl 4 /c 3 



spin ■ 



is proportional to the spin-down luminosity, L. 
Therefore, the evolution law, L 7 oc L^ in , is crucially 
governed by the evolution of h m as a function of the 
NS age, t. 

The evolution of h m is essentially controlled by the 
photon-photon pair production in the pulsar magneto- 
sphere. To analytically examine the pair production, 



we assume the static dipole magnetic field configura- 
tion for simplicity, and consider the plane on which 
both the rotational and magnetic axes reside. On 
this two-dimensional latitudinal plane, the last-open 
field line intersects the NS surface at magnetic co- 
latitudinal angle #" lax (measured from the magnetic 
dipole axis) that satisfies 



sin' 



; (0lc-") 



tJ7Lc/sin6>LC ' 



(5) 



where r* denotes the NS radius, #lc the angle (mea- 
sured from the rotation axis) of the point where the 
last-open field line becomes tangential to the light 
cylinder, and a the inclination angle of the dipole 
magnetic axis with respect to the rotation axis. A 
magnetic field line can be specified by the magnetic 
co-latitude (measured from the dipole axis), 6**, where 
it intersects the stellar surface. A magnetic field line 
does not close within the light cylinder (i.e., open to 
large distances) if < 8* < #" lax . Thus, the last-open 
field lines, 6** = 6>™ ax , corresponds to the lower bound- 
ary, which forms a surface in a three-dimensional mag- 
netosphere, of the outer gap. 

Let us assume that the gap upper boundary coin- 
cides with the magnetic field lines that are specified 
by 6** = (1 — /i m )T ax - Numerical examinations show 
that h m , indeed, changes as a function of the dis- 
tance along the field line and the magnetic azimuthal 
angle (measured counter-clockwise around the dipole 
axis). Nevertheless, except for young pulsars like the 
Crab pulsar, an assumption of a spatially constant 
h m gives a relatively good estimate. Thus, for an an- 
alytical purpose, we adopt a constant h m in this an- 
alytical examination. In this case, we can specify the 
middle-latitude field line by the magnetic co-latitude 
0* = (1 - h m /2)6? ax . Screening of E {{ due to the po- 
larization of the produced pairs, takes place mostly 
in the lower altitudes. It is, therefore, appropriate to 
evaluate the screening of £j| around the point (ro,0o) 
where the null-charge surface intersects the middle- 
latitude field line. 

An inwardly migrating electron (or an outwardly 
migrating positron) emits photons inwards (or out- 
wards), which propagate the typical distance l± (or 
I2) before escaping from the gap. Denoting the cross 
section of the inward (or outward) horizontal line from 
the point (ro,#o) and the upper boundary as (r\,Q\) 
(or as (r2,0 2 )), and noting r cos#o = ?"iCos#i = 
r2 cos 62 , we obtain 



l\ = ?*o cos 8q (tan 9q — tan 9\ ) , 



(G) 



h — To cos 6*o (tan 6*2 — tan# )- (7) 
Along the upper-boundary field line, we obtain 
sin 2 (6»i -a) sin 2 (6» 2 -a) sin 2 [(l - /i m )(9™ ax ] 



'2 



(8) 
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whereas along the middle-latitude field line, we obtain 
sin 2 (# -a) sin 2 [(l-/i m /2)^] 



ro 



(9) 



Combining these two equations, and noting 0™ ax <§; l ; 
we find that 8\ (< 8 ) and 8 2 (> 8 ) can be given by 
the solution 8 that satisfies 



cos 8 sin 2 (# — a) 



where 8q is given by 



1 - h m 
1 - h m /2 



cos 8o sin 2 (8$ — a) , 
(10) 



tan 8q — — (j> tan a + \J 9 tan 2 



(11) 



Thus, if we specify a, we can solve 8 = 8\ and 8 — 82 
as a function of h m by equation (1101) . Substituting 
these 8\ and 6*2 into equations © and (JT)), we obtain 
h and Z2, where ro depends on P. 

If hm <C 1, we can expand the left-hand side of equa- 
tion (fTU)) around 8 — 80, where 8 = 8\ for inward (or 
8 = 82 for outward) 7-rays to find 8 2 — 8 = 80 — 8\ oc 
>/h m . That is, the leading terms in the expansion van- 
ish and we obtain 1% = l 2 from the next-order terms, 
which are quadratic to 8 — 8$. Assuming Lx oc t~P , 
where f3 w 0.48 is appropriate for t < 10 4 years for 
a light-element-envelope NS and for t < 10 5 years for 
a heavy-element NS, we find h m oc F 5 / 6 w" 1/6 t /3/2 , 

v 



and hence L, k F 3/2 M 3/2^3/3/2_ Since the dipole 



radiation formula gives P oc /it 1 / 2 , we obtain L 7 oc 
/i o^3(/3-i/2)/2 K i -o.03_ ThuSj when the gap ig very 

thin, which is expected for a light-element-envelope 
NS, L 7 little evolves with the pulsar age, t. 

However, if h m > 0.2, say, the rapidly expand- 
ing magnetic flux tube gives asymmetric solution, 
I2 > h- That is, the third and higher order terms in 
the expansion contribute significantly compared to the 
quadratic terms. Thus, we must solve equation (fTOj) 
for 8 (= 6*i or 82) without assuming \6 — 9q\ <C 1, in 
general. 

Let us now consider the condition for a gap to be 
self-sustained. A single ingoing e~ or an outgoing e + 
emits 



or 



(JV 7 )! = eE\\h/{hv c ) 



(Ag 2 = eE\\l 2 /{hv c ) 



(12) 



(13) 



photons while running the typical distance or l 2l 
respectively. Such photons materialize as pairs with 
probability 



or 



T2 = I2F2CF2/C, 



(14) 



(15) 



where F\ and F 2 denote the X-ray flux inside and 
outside of (r ,# ), respectively; o~\ and U\ are the pair- 
production cross section for inward and outward pho- 
tons, respectively. Thus, a single e~ or e + cascades 
into 



pairs or into 



eE llFl 2 

ti G\ 

hv c c 



eE {] F 2 2 

1 2 02 



(16) 



(17) 



pairs within the gap. That is, a single inward- 
migrating e~ cascades into pairs with multiplicity 
(iV 7 )iTi. Such produced pairs are polarized by E». 
Each returning, outward-migrating e + cascades into 
pairs with multiplicity (A 7 )2T2 in outer magneto- 
sphere. As a result, a single inward e~ cascades even- 
tually into (iV~)iTi'(iV~)2T2 inward e~'s, which should 
become unity for the gap to be self-sustained. There- 
fore, in a stationary gap, the gap thickness h m is auto- 
matically regulated so that the gap closure condition, 



(7V 7 ) 1 r 1 • (7V 7 ) 2 r 2 - 1, 



(18) 



may be satisfied. 

Approximately speaking, a single, inward-migrating 
e~ emits (AT 7 )i ~ 10 4 curvature photons, a portion of 
which head-on collide the surface X-ray photons to 
materialize as pairs with probability n ~ 10 -3 within 
the gap. Thus, a single e~ cascades into (iV~)iri ~ 10 
pairs in the gap. Each produced e + return outwards to 
emit {N- l )2 ~ 10 5 photons, which materialize as pairs 
with probability t 2 ~ 10~ 6 by tail-on colliding with 
the surface X-rays. In another word, (N 7 )iTi ~ 10 
holds regardless of the nature of the pair produc- 
tion pro cess (e.g., either p hoton- photon or magnetic 
process [Takata et al.l l2010l |) in the lower altitudes, 
because it is determined by the pair-production ef- 
ficiency in the outer magnetosphere (-/V 7 ) 2 t 2 ~ 0.1, 
which is always due to photon-photon pair produc- 
tion. 

In general, (-/V 7 )i, n, (-^7)2, t 2 are expressed in 
terms of h m , P, /i, T, and a, where T denotes the 
NS surface temperature. Note that we can solve P = 
2ir/Vl as a function of the NS age, t, from the spin- 
down law. Thus, specifying a and the cooling curve, 
T = T(t), we can solve h m as a function of t from the 
gap closure condition, (A r 7 )iri(A r 7 ) 2 T 2 = 1. Note also 
that the spin-down law readily gives the spin-down 
luminosity, L sp ; n oc PP~ 3 , as a function of t, once 
P = P{t, a) is solved. On these grounds, L 7 can be 
related to L sp - ln with an intermediate parameter t, if 
we specify the cooling curve and the spin-down law. 

Substituting equations p^|) and ([!?)) into ([18"]). we 
obtain 



e£j| ^F 1 <t 1 F 2 (T2 
hv c c 



hh = 1, 



(19) 
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where 



( \ 2 />oo 
~ ) / 
r */ J» th , 



BAT) 
hv 



(20) 

with i = 1,2; v 1 denotes the 7-ray frequency, and 
B v (T) the Planck function. We have to integrate over 
the soft photon frequency v above the threshold en- 
ergy 



2\2 



hi/tb,- 



2(m c c 2 ) 
(1 - Hi)hv 1 



(21) 



where m c c 2 refers to the rest-mass energy of the elec- 
tron. The cosine of the collision angle fa becomes 
1 — fix = 1 — sin 8q for outward (or 1 — fj,2 = 1 + sin#o 
for inward) 7-rays. That is, collisions take place head- 
on (or tail-on) for inward (or outward) 7-rays. The 
total cross section is given by 



cr P = — cr T (l - v 2 ) 
16 



1 — V 



2v(2-v 2 ) 



where ctt denotes the Thomson cross section and 



(22) 



1- 



(m e c 2 ) 2 



1 — ^ hvhv-f 



(23) 



Pair production takes place when the 7-rays col- 
lide with the surface X-rays in the Wien regime, that 
is, at hv 3> kT. An accurate evaluation of cr 2 re- 
quires a careful treatment of the collision geometry, 
because the threshold energy, hv t h,2, strongly depends 
on the tiny collision angles. In the numerical method 
(next section), the pair-production absorption coeffi- 
cient is explicitly computed at each point in the three- 
dimensional pulsar magnetosphere. However, in this 
section, for analytical purpose, we simply adopt the 
empirical relation, 



(24) 



where e = 0.004, 0.01, and 0.038 for a = 45°, 60°, 
and 75°, respectively; 1 — fii ~ 2. The X-ray flux is 
evaluated at (ro,#o) such that 



F x 



1 



2.82fcT 47rr 2 ' 



(25) 



where Lx refers to the luminosity of photon radiation 
from the the cooling NS surface. For a smaller a, the 
point (^,#2) is located in the higher altitudes, where 
the magnetic field lines begin to collimate along the 
rotation axis, deviating from the static dipole configu- 
ration. Thus, the collision angles near the light cylin- 
der, and hence 02 decreases with decreasing a. The 
explicit value of e can be computed only numerically, 
solving the photon specific intensity from infrared to 



7-ray energies in the three-dimensional pulsar magne- 
tosphere. 

The last factor, Ziig, in the left-hand side of equa- 
tion (|19[) is given by 

hh = r Q 2 cos 2 o (tan0 o — tan0i)(tan0 2 — tan0 o ) (26) 



= (27) 



Thus, equation (|T9|) gives 
eE\\ L x /c 



hv c 2.82/fcT 
x cos 2 



'0 (tan 6*o — tan 8\) (tan #2 — tan#o) 



where the ro dependence vanishes. Substituting equa- 
tions (P) , ([!]) , © into (|2T| , we can solve h m as a func- 
tion of Lx/kT, P, and /x. 

To describe the evolution of P = P(t) = 2n/Q(t), 
we adopt in this paper 



- mn = c 



H 2 n 4 



(28) 



where C — (2/3) sin 2 a for a magnetic dipo le braking, 
while C = 1+sin 2 a for a force- free braking [Spitkovskil 
l2006j . Assuming a magnetic dipole braking, we obtain 



P 



■ 39.2ms/i 3 o^45 1/2 (t/10 3 years) 1/2 , 



(29) 



where ^ 30 = ^/(10 30 Gcm 3 ) and Z45 = 
I/(10 45 g cm 2 ). Thus, if we specify a cooling 
scenario, T = T(t), equation (f2"7| gives h m as a 
function of t. Note that the a dependence of the 
spin-down law is not essential for the present purpose; 
thus, C = 2/3 is simply adopted. Once /i m = h m (t) is 
obtained, equation ^ readily gives L 7 as a function 
of t, and hence of L sp i n . It is worth noting that 
the heated polar-cap emission is relatively weak 
compared to the cooling NS emission, except for 
millisecond or middle-aged pulsars. 

We adopt the minimal cooling scenario [Page et ail 
2004] Within the minimal cooling scenario, the cooling 
history of a NS substantially depends on the compo- 
sition of the envelope. W e adopt the cooling curves 
given in Pag e et aTl [2004] and consider the two ex- 
treme cases: light element and heavy element en- 
velopes. 

We present the solved h m for a light a n d a h eavy 
element envelope in figure 2 in iHirotanil [2013j . It 
follows that the gap becomes thinner for a light el- 
ement case than for the heavy element cases. This 
is because the more luminous photon field of a light 
element envelope leads to a copious pair production, 
which prevents the gap to expand in the trans-field 
direction. As a result, the predicted L 1 becomes less 
luminous for a light element envelope than a heavy 
one. 

In fi gure [Jl we present the analytical results of 
versus L S pin as the dotted (or dashed) curve for a light 
(or a heavy) element envelope. As the pulsar spins 
down, L 7 evolves leftwards. It is interesting to note 
that L 7 littl e evolves for a light element envelope, as 



explained in IHirotanil [2013 
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3. Numerical examination of outer-gap 
electrodynamics 



Let us develop the analytical examination and look 
deeper into a self-consistent solution by a numer- 
ical method. T o this end, we adopt the modern 
outer-gap model [Hirotanil 1201 laj and solve the set 
of Maxwell and Boltzmann equations self-consistently 
and compute E\\, distribution functions of e ± 's, and 
the photon specific intensity at each point in the three- 
dimensional pulsar magnetosphere. We consider not 
only the whole-surface, cooling NS emission but also 
the heated polar-cap emission as the photon source of 
photon-photon pair production in the numerical anal- 
ysis. The former emission component is given as a 
function of the pulsar age from the minimum cool- 
ing scenario, in the same manner as in the analytical 
examination, while the latter emission component is 
solved consistently with the energy flux of the e~'s 
falling on to the pulsar polar-cap surface. 

The numerical method is described in IHirotanil 
2013] in detail. We solve the set of partial and or- 
dinary differential equations under the boundary con- 
ditions that e^'s or 7-rays do not penetrate into the 
gap from outside. By this method, we can solve 
the acceleration electric field En, particle distribu- 
tion functions n±, and the photon specific intensity 
l v (from hv = 0.005 eV to 50 TeV), at each position 
in the three-dimensional magnetosphere of arbitrary 
rotation-powered pulsars, if we specify P, ji, a, and 
kT. We adopt the minimum cooling scenario in the 
same manner as in § [2] 

In figure [TJ we plot the result of L 7 as a function of 
L S pin as the dash-dotted (or solid) curve for a light (or 
a heavy) element envelope, where /X30 = 3.2 is adopted 
in the same manner as in the analytical examination. 
It follows that these numerical solutions are consis- 
tent with the analytical ones, and that L 7 decreases 
slowly until 10 4 ' 5 years. The physical reason why L 1 
increases with decreasing L sp j n at t > 10 4 years for 
a light element envelope, is the same as described 
at the end of § [2] A realistic NS will have an en- 
velope composition between the two extreme cases, 
light and heavy elements. Thus, the actual L 7 's will 
distribute between the red solid (or dashed) and the 
blue dash-dotted (or dotted) curves. However, after 
.L-y approaches £ S pi n (t hin dashed straight line; see 
I Wang fc Hirotanil 201 lj for the death line argument), 
the outer gap survives only along the limited magnetic 
field lines in the trailing side of the rotating magne- 
tosphere because of a less efficient pair production; as 



a result, L 7 rapidly decreases with decreasing L S pin- 
For a smaller a, even for a light element envelope, 
L 1 monotonically decreases as the dash-dot-dot-dot 
curve shows, because the gap is located in the higher 
altitudes, and because the less efficient pair produc- 
tion there prevents the produced electric current to 



increase with decreasing age around t ~ 10 years 



4. Discussion 

To sum up, a light element envelope approximately 
corresponds to the lower bound of the (observationally 
inferred) gamma-ray luminosity of rotation-powered 
pulsars, whereas a heavy element one to the upper 
bound. The scatter of the intrinsic gamma-ray lumi- 
nosity is physically determined by the magnetic incli- 
nation angle, a, and the envelope composition. The 
cutoff energy of the primary curvature emission is kept 
below several GeV even for young pulsars, because the 
gap trans-field thickness, and hence the acceleration 
electric field, is suppressed by the polarization of the 
produced pairs in the lower altitudes. 

To convert the observed 7-ray flux into luminosity, 
Lj, one has conventionally assumed /n = 1. For ex- 
ample, the error bars of the observational data points 
in figure [TJ do not contain any uncertainties incurred 
by fn- Nevertheless, if a and £ can be constrained, 
we can estimate L 7 more accurately, by applying the 
present quantitative outer-gap calculations. It is note- 
worthy that L 7 's given in figure [TJ little depend on the 
NS magnetic moment, /1. This is particularly true for 
a light element case, which has h m <C 1, by the reason 
described after equation (jTTJ) . What is more, with an 
additional determination of d (e.g., by parallax obser- 
vations), we can infer the composition of individual 
NS envelopes, by using the constraine d flux correc- 
tion factor, /n (fig. 6 in IHirotanil |2013j |). We hope to 
address such a question as the determination of a and 
C, and hence fn, for individual pulsars, by making an 
'atlas' of the pulse profiles and phase-resolved spectra 
that are solved from the basic equations in a wide pa- 
rameter space of P, fi, T, a, and £, and by comparing 
the atlas with the observations. 
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Pulsar spin-down luminosity, ^ S pj n (erg s ) 

Figure 1: Evolution of the outer-gap luminosity as a function of the neutron-star spin-down luminosity. Analytical 
results are plotted as dotted and dashed curves, while numerical ones as dash-dotted, solid, and dash-dot-dot-dot ones. 
A light or a heavy element envelope is assumed, as indicated in the box. For comparison, the case of a = 45° is also 
depicted. The green filled circles designates the normal gamma-ray pulsars, while the blue filled squares do those 
detected by the gamma-ray blind search technique. 
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